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Abstract. We present the first results of the analysis of data collected during the 1998-99 observational campaign at 
the 1.3 meter McGraw-Hill Telescope, towards the Andromeda galaxy (M31), aimed to the detection of gravitational 
^ ^, microlensing effects as a probe of the presence of dark matter in our and in M31 halo. The analysis is performed using 
the pixel lensing technique, which consists in the study of flux variations of unresolved sources and has been proposed 
and implemented by the AGAPE collaboration. We carry out a shape analysis by demanding that the detected flux 
variations be achromatic and compatible with a Paczyiiski light curve. We apply the Durbin- Watson hypothesis test 
C3 ' to the residuals. Furthermore, we consider the background of variables sources. Finally five candidate microlensing 
events emerge from our selection. Comparing with the predictions of a Monte Carlo simulation, assuming a standard 
spherical model for the M31 and Galactic haloes, and typical values for the MACHO mass, we find that our events 
are only marginally consistent with the distribution of observable parameters predicted by the simulation. 
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1. Introduction 

In the last decade much attention has been focused on 
the possibihty that a sizable fraction of galactic dark 
matter consist of MACHOs (Massive Astrophysical Com- 
pact Halo Object). Since 1992, the MACHO (Alcock et 
al. 1993) and EROS (Aubourg et al. 1993) collabora- 
tions have looked towards the Large and Small Magellanic 
Clouds (LMC and SMC) in order to detect MACHOs us- 
ing gravitational microlensing. This technique, originally 
proposed by Paczyhski (1986), analyses the luminosity 
variation of resolved source stars, due to the passage of 
MACHOs close to the line of sight between the source 
and the observer. 

The MACHO collaboration (Alcock et al. 2000) dis- 
covered 13 - 17 microlensing events towards the LMC. 
Assuming that all events are due to MACHOs in the halo, 
about 20% of the halo dark matter resides in form of com- 
pact objects with a mass in the range 0.15 — 0.9Mq. The 
EROS collaboration (Lasserre et al. 2000) observed 6 mi- 
crolensing events, 5 in the direction of the LMC and 1 
in the direction of the SMC. These observations place an 
upper limit on the halo dark matter fraction in the form 
of MACHOs. In particular, they exclude, at the 95% con- 
fidence level, that more than 40% of a standard halo is 
composed of objects in the range IO^^Mq — IMq. Note 
that the results of the two collaborations are consistent 
with a 20% halo dark matter fraction of objects ~ OAMq. 

The OGLE collaboration (Udalski et al. 1993) origi- 
nally searched for microlensing events only towards the 
Galactic bulge, but has now also extended its search to 
the LMC and SMC. 

A natural extension of the microlensing observational 
technique consists in observing dense stellar fields even if 
single stars cannot be resolved, as in the case of the M31 
galaxy. For this purpose, the pixel lensing technique has 
been proposed (Baillon et al. 1993) and then implemented 
by the AGAPE collaboration (Ansari et al. 1997). Another 
technique, based on image subtraction, has been devel- 
oped by the VATT-Columbia collaboration (Crotts 1992; 
Tomaney & Crotts 1994), and is used also in the WeCAPP 
project (Riffeser et al. 2001). The monitoring of M31 has 
the advantage that the Galactic halo can be probed along 
a line of sight different from those towards the LMC and 
SMC. Furthermore, the observation of an external galaxy 
allows one to study its halo globally, which, in the case 
of M31, has a particular signature due to the tilted disk. 
Accordingly, the expected optical depth for microlensing 
varies from the near to the far side of the M31 disk (Crotts 
1992; Jetzer 1994). 

The efficiency of the pixel lensing method to detect 
luminosity variations has been tested by the AGAPE col- 
laboration on data taken at the 2 meter Bernard Lyot 
Telescope in two bandpasses {B and i?), covering 6 fields 
of 4.5' X 4.5' each around the center of M31, during 3 years 
of observations (1994-1996). A possible microlensing can- 



didate has been observed and further characterized by us- 
ing information from an archival Hubble Space Telescope 
WFPC2 image (Ansari et al. 1999). An important con- 
clusion of this analysis is that it is crucial to collect data 
in two bandpasses over a long duration with regular sam- 
pling. Very recently, the POINT-AGAPE collaboration 
(Auriere et al. 2001) announced the discovery of a short 
timescale candidate event towards M31. Additional mi- 
crolensing candidates towards the same target have been 
reported by the VATT-Columbia collaboration (Crotts et 
al. 2000). 

In this paper, we present results for the 1998-1999 
campaign of observations at the 1.3 meter McGraw-Hill 
Telescope, MDM Observatory, Kitt Peak, towards the An- 
dromeda Galaxy. In § H, we briefly outline the pixel lensing 
technique. § § is devoted to the description of the obser- 
vational campaign and the experimental setup. In § || we 
discuss the data reduction procedure (Calchi Novati 2000) 
in some detail, in particular the approach used to elimi- 
nate instabilities caused by the seeing and to evaluate the 
errors. In § || we present our selection pipeline (Calchi 
Novati 2000): bump detection (§ |3.l[ ), shape analysis (§ 
|S.2[ ) and color and timescale selection (§ 5.3). We select 
a sample of 5 light curves that we retain as microlensing 
candidate events and whose characteristics are given in § 
.41 In § 5.5 we show the light curve of a nova located 



inside our field of observation: the discovery of variable 
sources is a natural byproduct of the microlensing search. 
In § ^ we conclude with a comparison of the outcome of 
our selection with the prediction of a Monte Carlo simu- 
lation. 

2. Pixel Lensing 

Pixel lensing is an efficient tool for searching microlens- 
ing events when the sources can not be resolved. In this 
case, the light collected by each pixel is emitted by a huge 
number of stars. Although, in principle, all stars in the 
pixel field are possible sources, one can only detect lens- 
ing events due either to bright enough stars or to high 
amplification events. Typical sources are red giant stars 
with Mj e [—3.5,0]. We estimate that there are about 
100 sources per square arc second that fit these require- 
ments. The main drawback of the method is that usually 
we have no direct knowledge of the flux of the unamplifled 
source. 

For images obtained by the observation of dense stel- 
lar fields, the flux collected by a pixel is the sum of the 
fluxes emitted by single stars, which all contribute to the 
background. If one of these stars is lensed, its fiux varies 
accordingly. Whenever this variation is large enough, it 
will be distinguishable from the background produced by 
the other stars. Denoting by (j) the amplified flux detected 
by the pixel and by (ptkg the background flux, the fiux 
variation is given by 



- ^bkg = r [Ait) - 1] 



(2.1) 
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Fig. 1. M31 with MDM and Agape observation fields (courtesy of A. Crotts). White stars and dots give, respectively, 
the position of the five microlensing candidate events (labelled as in Table ^ and where candidates 1 and 2 appear to 
be superimposed) and of the nova. 



where A{t) represents the amplification as a function of 
time. 



(f'bkg = 4>* + 4>res , 



(2.2) 



(jj* is the flux of the star before lensing, and (j)res is that 
given by the other stars. In the point-like and uniform- 
motion approximations, A is related to the lensing param- 
eters by: 



A 



2 _ (t- ^o)^ , „,2 



^2 



(2.3) 
(2.4) 



where is the Einstein time, uq the impact parameter 
in units of the Einstein radius i?g , and to the time of the 
maximum amplification. The Einstein radius is 



Re 



AGM DoiDi, 



(2.5) 



where M is the mass of the lens, Doi, Dos and Dis are the 
observer-lens, observer-source, and lens-source distances, 
respectively. 

Two characteristic features of a microlensing event are 
achromaticity and the uniqueness of its luminosity bump, 
although differential amplification of extended sources can 
give rise to a chromatic, but still symmetric, lensing light 
curve (Han et al. 2000). 
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Fig. 2. Time sampling for the observations in the "Target" 
field. 



Most of the observations (about 80%) are concentrated 
in the first three months, so that, unfortunately, the time 
distribution of the data is not optimized for the study of 
microlensing effects. Thus, the given time distribution al- 
lows us to select events that take place almost exclusively 
during the first three months of observation. Furthermore, 
the time coverage of about 14 months is still not long 
enough to test conclusively the bump uniqueness require- 
ment for a microlensing event. Mainly for this reason, we 
will speak in this paper only of candidate microlensing 
events. 

Taking into account the transmission efficiency of the 
filters and the catalogued magnitudes Rc and Ic (Cousins 
colour system), for a sample of 23 reference secondaries 
identified in the Target field (Magnier et al. 1993), we 
derive the following photometric calibrations 

Rc = ruR - 0.13 • {wiR - mi) + 22.54 , (3.6) 

/c = m/ - 0.02 -(mfl- TO/) + 22.21, (3.7) 

where niR^-^ = —2.5 • log(0^^^.j) and 0^^/^ is the flux of 
the source in ADU/s measured in the R and / filters re- 
spectively. The estimated error is ~ 0.1 mag, both for R 
and /. 



3. Observations and experimental setup 

The data analysed in this paper have been collected on 
the 1.3 meter McGraw-Hill Telescope, at the MDM obser- 
vatory, Kitt Peak (USA)Q Two fields have been observed, 
which lay on the two sides of the galactic bulge (see Fig. 
|l|) and have been chosen in order to be able to study the 
expected gradient in the optical depth. The two fields are 
almost parallel to the major axis of M31; their centers are 
located at a = 00h43m24s, S = 41° 12' 10" (J2000) (named 
"Target"), and a = 00h42ml4s, S = 41° 24' 20" (J2000) 
(named "Control"). The data acquired in the Target field 
are analysed here. 

Fig. |l| shows the location of the fields and for compar- 
ison also the smaller AGAPE field. The observations were 
taken with a CCD camera of 2048 x 2048 pixels with 0"5 
and therefore a total field size of 17' x 17'. 

In order to test for achromaticity, images have been 
taken in two bands, a wide R and a near--standard /. The 
exposure time is 6 minutes for R, 5 minutes for /. The 
observations started in the fall of 1998 and are still under- 
way. 

Here we analyse the data taken in the period from the 
beginning of October 1998 to the end of December 1999. 
In Fig. 1^ we give the time sampling of the measurements 
(number of nights and images). For each night we have 
measurements in both R and /. On average, there are 
twice as many R images as / images. In R band we have 
~ 800 images distributed along 42 nights of observation. 

^ Data shared with the Columbia- Vatt collaboration. 



4. Data reduction 

During each night of observation about 20 images are 
taken in R and 12 in /. In principle, this allows two possi- 
ble strategies for the analysis. We can study fiux variations 
on light curves built either by a point obtained from each 
image, or by a point obtained by averaging over many 
images. In the first case, we are potentially sensitive to 
very short time variations. However, this sensitivity is un- 
dermined by the low signal-to- noise ratio (S/N). In the 
second case, the S/N is increased by the square root of 
the number of images we combine. 

Results of the analysis of light curves built with one 
point per image will be discussed in a future paper. Here 
we concentrate on the analysis of light curves obtained af- 
ter combining all the images taken in the same night, using 
a simple averaging procedure performed on geometrically 
aligned images (see below). 

Data reduction is carried out as follows. After the usual 
corrections for instrumental effects, debiasing and flat- 
fielding, we normalize all the images to a common refer- 
ence to cope with variations induced by the observational 
conditions which arc different from image to image (so 
that we get global stability conditions on each image with 
respect to a given one). We can distinguish three separate 
effects: the geometric offset of each image with respect to 
the others, the difference in photometric conditions of the 
sky and seeing effects. 

By means of geometrical alignment we obtain that each 
pixel, on all the images, is directed towards the same por- 
tion of M31. We take advantage of the fact that the mean 
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seeing disk is much larger than the pixel size. We follow 
Ansari et al. (1997), and get a precision better than 0."1. 

Following the methods developed by the AGAPE col- 
laboration (Ansari et al. 1997), we then bring all the im- 
ages to the same photometric conditions in such a way 
that the images are globally normalized to a common ref- 
erence. The procedure is based on the hypothesis that a 
linear relation exists between "true" and measured flux. 
It then follows that 



I pixel / ■ -\ 
Kef i^'j) 



Veurr J ) ' "c 



(4.8) 



where </)^gj'^'(z, j) and 0^™^'(z,j) are the fluxes in pixel 
for the reference and the current image, and Qcurr 
and bcurr are two correction coefficients, which are the 
same for all the pixels of the image, that take into account 
the effects of variable atmospheric absorption and variable 
sky respectively. 

The seeing effect gives rise to a spread of the received 
signal. In our data the seeing varies from 1"3 up to 
~ 2."2. Consequently we observe fake fluctuations on light 
curves obtained after photometrical alignment. In order 
to cope with this effect, and thus to get reasonably stable 
light curves, we follow a two steps procedure (for further 
details see Ansari et al. (1997) and Le Du (2000)). We 
begin by substituting the flux of each pixel with the flux of 
the corresponding superpixel, defined as the flux received 
on a square of m x m pixels around the central one. The 
value m should be large enough to cover the typical seeing 
disk, but not too large to avoid an excessive dilution of the 
signal. Given the mean seeing value and the angular size 
of the pixel, we choose m = 5. This corresponds to 2."5, 
compared to the average value for the seeing of ~ l."7 for 
both R and / images. In this way we get a substantial gain 
in stability since elementary pixels are strongly affected by 
seeing fluctuations. 

Denoting by $(«, j) the flux in a superpixel, we have 



i+n 



j+n 

E 

k—i—n l—j — n 



(4.9) 



where n = (to — l)/2 and to is the superpixel size. 

Instead of trying to evaluate the point spread func- 
tion of the image, as a second step we apply an empirical 
stabilization of the difference between the flux measured 
on the image and that of the median image, obtained by 
removing small scale variations with a median filter on 
a very large window of 31 x 31 pixels (in this way we 
get an image whose signal is independent from the seeing 
value). The stabilization is then based on the observed 
linear correlation, for each superpixel, between these dif- 
ferences measured on the current image (after photomet- 
rical alignment) and the reference image. Denoting with 
$,.e/(i, j), ^curr,a{i,j) and $"'='*(«, j) the value of the flux 
in a superpixel for the reference, the current (pho- 

tometrically aligned) and the median images respectively. 



1.3 1.4 1.5 1.6 1.7 1.8 1.9 2.0 



seeing (arcsec) 
Fig. 3. The value of the correction factor 

c^curr in the 



relation (4.10) for each composed image as a function of 
the seeing. 




Fig. 4. Plot showing the linear correlation be- 
tween the quantities [^curr,aii, j) — ^™urr.a(h j)] and 



$,.e/(*,j)-<i>™//(i,j) 

line. 



The dashed line is the y 



we have the empirical relation 

'^curr,a{h j) curr.ai^ ^ j) ~ 



(4.10) 



The slope acum calculated with a minimization proce- 
dure, shows a clear correlation with the seeing (Fig. ||). 
This variation is expected because the flux of a star that 
enters a superpixel changes with the seeing. 

In the example shown in Fig. ^jthe seeing of the current 
image is greater than that of the reference image, and we 
find acurr < 1- We note that in this case, the flux of 
superpixels for which ^curr,a < ^™urr a are corrected to 
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Fig. 5. The same light curve before (top) and after see- 
ing correction. In both cases the error bar shows just the 
photon noise. 



a lower value (points in the bottom left corner), while, if 
^curr,a > '^Turr ^hey are corrected up (points in the up 
right corner). If the seeing of the current image is smaller 
than that of the reference image {acurr > 1) the situation 
is reversed. 

While stressing its empirical character, we note that 
this approach is rapid and efficient. In Fig. ^ we show the 
effect of the correction on a given light curve. 

To construct a corrected current pixel flux as close as 
possible to the reference flux, $re/, we replace $re/ in 
(4.10) by the corrected current flux ^curr and solve for 
the latter: 



(4.11) 



That is, the corrected flux for the current image is given 
by the sum of the median of the reference image, and of 
the weighted deviation from its median, i. e. by its char- 
acteristic small spatial scale variations, depending on the 
relative absorption and on the seeing conditions^. 

In order to minimize the deviations from the median, 
we choose as a reference image (one for each filter), an 
image characterized by a seeing value equal to the average 
value over the period of observations. The seeing fraction 
of a source in a superpixel is then ~ 0.87. 

Another crucial point is the evaluation of the error 
to be associated with the received flux. In order to give 
a more appropriate, though empirical, error estimate, we 
renormalize the photon noise (Tstat by introducing a correc- 
tion factor rjcurr, that depends on the image, to include all 
the systematic effects over which we have less control. For 



W e have verified that it is possible to exploit the relation 
(4.10) with ^curr.a replaced by $curr (i.e. using non photomet- 



rically aligned images) in order to perform in one single step 
the photometrical alignment and the seeing correction. In the 
two cases we get the same final results. 



1.3 1.4 1.5 1.6 1.7 1.8 1.9 



: 2.1 2.2 2.3 

seeing (arcsec) 



1.3 1.4 1.5 1.6 1.7 1.8 1.9 



: 2.1 2.2 2.3 

seeing (arcsec) 



seeing (arcsec) 



Fig. 6. The dispersion for the distribution ( 4.12 ) calcu- 
lated for each composed image as a function of the seeing 

before (^(ycurrjbefore 

) and after {acurr, after) seeing correc- 
tion (taking all the points of the image) and for the sample 
of points selected according to the condition (4.15). 



a discussion on the relation between the photon noise and 
others different systematic effects such as surface bright- 
ness fluctuations see Gould (1996). 

The evaluation of rjcurr is based on the study of the 
dispersion of the distribution of the normalized difference, 
superpixel by superpixel, between the current and the ref- 
erence image 



(4.12) 



^currjij) - ^refjij) 
a Stat curr 

where the denominator is the statistical error on the dif- 
ference. 

This distribution is expected to have zero mean (which 
follows from the geometrical and photometrical align- 
ment) and dispersion one (which would indicate that the 
photon noise alone gives the right evaluation of the error). 
We do find a null mean but the dispersion is greater than 
one and depends on the seeing value. We note, however, 
that this effect is greatly reduced by the seeing correction 
(see Fig. ^). For the reference image we assume that the 
estimated error (Test [$re/] coincides with the statistical er- 
ror given by the photon noise. 



(4.13) 



because it does not suffer from additional noise due to 
geometric, photometric and seeing transformations, as do 
the other images. On the other hand, this means that we 
evaluate the estimated error with respect to that of the 
reference image. 

The correction factor rjcurr is then equal to the disper- 



sion of the distribution ( 4.12 ) calculated for a sample of 
points, properly selected according to the criterion that 
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they belong to "stable" light curves, in order to exclude 
light curves which show real stellar flux variations. If the 
subset is small enough and homogeneous, the correction 
factor appears to be, as expected, almost independent of 
the seeing value, and has an average value rjcurr ~ 1-4 
(Fig. |[ bottom). The estimated error of the (i, j) super- 
pixel is obtained from the relation 



2 

^est 



The "stable" pixels are selected by imposing the con- 
dition 



(4.15) 



where $(i,j) is the average flux along the light curve, 
^bkgihj) is the baseline level and cTstati^iiJ) - ^bkgih j)] 
is the error associated with the evaluation of the averages. 
The baseline is defined as the minimum value taken by the 
average flux calculated along q consecutive points on the 
light curve. We take g = 6, in order to avoid underestima- 
tion of the baseline due to fluctuations. Choosing e ^ 1.5, 
the selected points are ^ 3% of the total. 

5. Candidates selection 

A microlensing event is characterized by specific features 
that distinguish it from other, much more common, types 
of luminosity variability, the main background to our 
search. In particular for a microlensing event the bump 

~ does not repeat; 

— follows the (symmetric) Paczynski shape; 

— is achromatic. 

Variable stars usually show multiple flux variations 
and have an asymmetric chromatic shape. Moreover, dif- 
ferent classes of variable stars are characterized by specific 
features, such as timescale variation and color, that can 
be used to distinguish them from real microlensing events. 

In the following we devise selection techniques that 
make use of these characteristics while taking account of 
the specific features of our data set. 



that the points follow a gaussian distribution around the 
baseline, we use the estimator L, the likelihood function, 
to measure the statistical significance of a bump. We want 
to give more weight to points that arc unlikely to be found, 
so that we define L as |^ 

i = - In (njeb„„pF($|$ > $j)) given ^tkg, , (5.16) 
where 



(4.14) p($|$ > 



1 



exp 



{^-'^bkgf 



.(5.17) 



L is then a growing function of the unlikelihood that a 
given variation is the product of random noise. This esti- 
mator is different from the usual definition leading to a ^ 
with n points, and has the advantage to give weight only 
to the positive deviations above threshold, which arc the 
ones of interest. 

For each light curve we denote by Li and L2 the two 
largest deviations, respectively. We fix a threshold Ltresh^ 
and we require Li > Ltresh to distinguish real variations 
from noise. Moreover, we fix an upper limit to the ratio 
L2/L1 to exclude light curves with more than one signif- 
icant variation. The shape analysis is then carried out on 
the superpixels that have the highest values of L in their 
immediate neighborhood since we find a cluster of pixels 
associated with each physical variation. This method suf- 
fers from a possible bias introduced by an underestimation 
of the baseline level (which we further analyse in the next 
section). 

We have carried out a complete analysis selecting the 
pixels with the following criteria: 

— exclusion of resolved stars; 

— Ltresh — 100; 

— L2/L1 < 0.1. 

This selection is made only on R images in order to reduce 
contamination by variable sources. In this way we take 
advantage of the fact that most luminous variables (to 
which we are anyway sensitive) show stronger variations 
in the / than R band. 

By using these peak detection criteria, the number of 
superpixels is reduced from ~ 4 • 10^ to ~ 5 • 10^. 



5.1. Bump detection 

As a first step we select light curves showing a single flux 
variation. We begin by evaluating a baseline, i.e. the back- 
ground flux ($bfcg) along each light curve, as defined in § 

i 

Once the baseline level has been fixed, we look for a sig- 
nificant bump on the light curve. This is identified when- 
ever at least 3 consecutive points exceed the baseline by 
3 (T. The variation is considered to be over when 2 consec- 
utive points fall below the 3 a level. Under the hypothesis 



5.2. (Achromatic) shape analysis 

As a second step we determine whether the selected flux 
variation is compatible with a microlensing event. 

The light curve of a microlensing event with amplifica- 
tion A{t) due to a source star with unlcnscd flux 0* (now 
to be evaluated in a superpixel) is 



Ht) = '^bkg + iA{t) - 1) r 



(5.18) 



^ In order to simplify the notation, hereafter we write $ fo r 
the corrected superpixel flux $curr as given in equation 4.11. 
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where represents the flux collected in the superpixel 
associated with a single pixel, as defined before, and A{t) 
is given by (2.3). 

Actually, one can not directly and easily measure 0*, 
the unamplified flux of the unresolved source star. Only 
a combination of the 5 parameters that characterize the 
light curve can be measured in a straightforward manner: 

^ ^bkg, the background level (which include the flux of 

the unamplified source); 
— to, the time of maximum amplification; 



t 



1/2 = ti/2 (^E 



'-J-min), the time width of the bump at 

half-maximum; 

— A^rnax = A^.^^^ {(/>*, Umin) = (f>* (^maa; - 1), the CX- 

cess of the flux with respect to the background at max- 
imum. 

Whenever the amplification is high enough, one can 
approximate A{t) ~ 1/u {t) and — 2^/3 Umin- It is 
then possible to rewrite the expression ( ^.1^ ) in terms of 
these 4 parameters, and a degeneracy arises among the pa- 
rameters of the amplification uo and tg, and the unknown 
flux of the unamplified source, 0* ("degenerate" Paczyhski 
curve Gould (1996)). Because of this degeneracy it is in 
general difficult to get, without extra-information, a reli- 
able insight into all the parameters that characterize the 
light curve, the Einstein time in particular. For this reason 
we can extract only the 4 aforementioned parameters even 
though we carry out a non linear fit with the complete 5 
parameters ("non-degenerate") Paczyhski curve. 

We now refine the selection based on the likelihood 
estimator in order to remove unwanted light curves with 
low S/N and for which the available data do not allow us 
to well characterize the bump. To this end we perform a 
Paczyhski 5-parameters fit and we study 

— the signal to noise ratio for the R flux variation; 

— the sampling of the data on the bump. 



We deflne the S/N estimator as 



\const 



Xii/n.d.f 



(5.19) 



where Xconst is the with respect to a constant flux and 
X^i is the with respect to the Paczyhski flt. 

The ratio Q is actually correlated with the likelihood 
estimator L we used in the previous step. In parallel with 
the cut Li > 100 we then keep only light curves with 
Q > 100. 

We do not ask for the / bump to be significant. 

The second point concerns the necessity to well char- 
acterize the bump shape in order to recognize it as a mi- 
crolensing event in the presence of highly irregular time 
sampling of data (see Fig. ||). For this purpose we require 
at least 4 points on both sides of the maximum, and at 
least 2 points inside the interval to ± ti/2/2. 

After this selection we are left with 1356 flux varia- 
tions. 



From now on we work with the data in both colors (J? 
and /) and we carry out a shape analysis of the light curve 
based on a two steps procedure as follows: 

— selection criterion; 

— Durbin- Watson test on residuals; 

which we now discuss in some detail. 

The first point is taken into account by performing 
the non-degenerate Paczyhski fit in both colors simulta- 
neously, so that we check also for achromaticity of the 
selected luminosity variations. In particular, we require 
that the three geometrical parameters that characterize 
the amplification (tg, uq and to) be the same in both col- 
ors. We get, therefore, a 7 parameters least non linear 
fit: 



2 Nj 

EE 

j=i i=i 



[$j(t,) - $We/(t,|$^fcg,(/'*'^to,tE'"o) 



.(5. 



To retain a light curve as a candidate microlensing 
event, we require that the reduced x'^ 



X 



N-7 



= f< 1.5 



(5.21) 



where N = Nn + Nj is the total number of points in / 
and R. 

The application of the criterion test reduces the 
sample of light curves from 1356 to 27. 

As a further step we apply the Durbin- Watson 
(Durbin & Watson 1951) test to the residuals with re- 
spect to the 7-parameters non-degenerate Paczyhski fit. 
With the DW test we check the null hypothesis that the 
residuals are not timely-correlated by studying possible 
correlation effects between each residual and the next one 
against type I error (i.e. against the error to reject the 
null hypothesis although it is correct, e.g. Babu & Feigel- 
son (1996)). We require a significance level of 10%. As 
time plays a fundamental role in the DW test, we perform 
this test on each color separately. 

We call dwn and dwi the coefficients for the Durbin- 
Watson test on the full data set. In order to retain a light 
curve we require 1.54 < dwfj(^j^ < 2.46, appropriate for 
40-42 points along the light curve. 

This statistical analysis reduces our sample of light 
curves from 27 to 11. 

It is worthwhile to note that some light curves, show- 
ing a real microlensing event superimposed on a signal 
due to some nearby variable source, and passing the pre- 
vious selection criteria, could be excluded by the DW test, 
sensitive to timely correlated residuals. 

In order to test our efficiency with respect to the in- 
troduction of the DW test we have done a study on "flat" 
light curves performing a "constant flux" fit, and selecting 
light curves requiring a reduced x^ < 1.5. By applying the 
DW test we then reject about 50% of these light curves. 
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i.e., more than the 10% we could expect if we had just | 
random fluctuations. In the discussion of the Monte Carlo « 
simulation we duly take into account this effect. 

•J 16D0 

I5Q0 

5.3. Color and timescale selection 



By far, the most efficient way to get rid of multiple flux 
variations due to variable stars is to acquire data that are 
distributed regularly for a sufficiently long period of time. 
Unfortunately, at present, the data cover with regularity 
only the flrst three months of observation, and the total 
baseline is less than 2 years long. 

For this reason, in addition to the analytical treat- 
ment that looks for the compatibility with an achromatic 
Paczyhski light curve (efficient, for instance, to reject 
nova-like events), we introduce another criterion based on 
the study of some physical characteristics of the selected 
flux variations. 

In particular we note that long period red variable 
stars (such as Miras) could not be completely excluded by 
the selection procedure apphed so far. By contrast, short 
period variable stars are eliminated thanks to the cut on 
the second bump of the likelihood function. A prelimi- 
nary analysis of the period, the color and the light curves 
of long period red variable stars, taken from de Laverny et 
al. (1997), lead us, with a rather conservative approach, to 
exclude those light curves that present at the same time 
a duration ti/2 > 40 days and a color {R — I)c > 1.0. 
A more detailed analysis aimed at a better estimation of 
that background noise is currently underway. 

We are aware that this last selection criterion could 
eliminate some real microlensing events. The probability 
that this might happen is however low because the mi- 
crolensing timescales are expected to be uncorrelated with 
the source color. In fact, a combination of a large ti/2 and 
color {R — I)c > 1.0 is extremely unlikely for microlensing 
events, but quite common for red variables. 

This last criterion further reduces the number of can- 
didates from 11 down to 5. 

5.4. Results of microlensing search 

We now summarize (Table |^) the different steps in the 
selection, give the number of surviving pixels after the 
application of the indicated criteria and the details of our 
set of microlensing candidate events. 

We take these 5 light curves, whose characteristics we 
are now going to discuss, as our final selection of mi- 
crolensing candidate events. 

In the following table (Table |^) we give their posi- 
tion, the estimated ti/2 in days, the time of maximum 
amplification (J-2449624.5), the magnitude at maxi- 



days 



< 705O 
2000 
1950 
19Q0 
1850 
IBOO 
1750 



davs 

Fig. 8. The R and / light curve of the nova event in 
a=00h43ml.7s,(5 = 41° 15' 37" (J2000). 

mum0 Rmax and the color {R — I)c. We then give the 
values of the fit: the reduced and the values of the 
Durbin-Watson dwfi(j) coefficients. 

The corresponding light curves are shown in Fig. 0. 

5.5. Variable stars 

Our data contain many more varying light curves that 
are due not to microlensing events but to other variable 
sources. The study of these variable stars is an interesting 
task in itself. Clearly, pixel lensing is well suited for this 
research. We give here (see Fig. |8|) only the light curve 
of an event characterized by a very strong and chromatic 
flux variation, that has already been considered as due to 
a nova (Modjaz & Li 1999) and whose light curve is also 
shown in Riffeser et al. (2001). We note that our points, 
even if with a poorer sampling in the descent, are in good 
agreement with those of this collaboration. This nova is 
found in the region of the bulge of M31 and we localize 
it in a=00h43ml.7s,(5 = 41° 15' 37" (J2000). We evaluate 
the magnitude at maximum amplification as Rc — 17.0 
and Ic = 16.8. 

6. Discussion of microlensing search 

In order to gain an insight into the results obtained, we 
compare our 5 candidate microlensing events with the pre- 
diction of a Monte Carlo simulation which takes into ac- 

^ Evaluated starting from the excess of flux with respect to 
the background A^rnax- 
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Fig. 7. Light curves of the 5 candidate microlensing events. On the y axis, A$ = $ — ^bkg- On the x axis, the origin 
of time is in J-2449624.5. The dashed hne represent the result of the 7-parameters Paczynski fit. 
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criterion 


% pixel excluded 


pixel left 


exclusion of resolved stars 


~ 10% 


~ 3.6 ■ 10" 


mono bump likelihood analysis 


~ 99.8% 


5269 


signal to noise ratio (Q > 100) 


~ 69% 


1650 


sampling of the data on the bump 


~ 18% 


1356 


X" < 1-5 


~ 98% 


27 


1.54 < dwR(i) < 2.46 


~ 59% 


11 


< 40 d or i? - J < 1.0 


~ 55% 


5 



Table 1. Summary of selection criteria. 



id 


q(J2000) 


5(J2000) 


ii/2 (d) 


to (J-2449624.5) 


Rmax 


R-I 




dwR 


dwi 


1 


00h43m27.4s 


41° 13' 11" 


32 ±6 


1506 ± 1 


22.6 ± 0.1 


0.3 ±0.2 


1.25 


1.78 


1.65 


2 


00h43m26.5s 


41° 13' 16" 


22 ± 7 


1508 ± 1 


22.7 ±0.1 


0.2 ± 0.2 


1.37 


1.57 


1.65 


3 


00h43m39.9s 


41° 18' 41" 


39 ±9 


1505 ± 1 


22.2 ± 0.1 


0.8 ±0.2 


1.48 


1.97 


1.67 


4 


00h42m39.3s 


41° 6' 53" 


15 ± 1 


1470 ± 1 


22.3 ± 0.1 


0.5 ±0.2 


1.42 


1.68 


1.95 


5 


00h42m39.1s 


41° 11' 26" 


25 ±3 


1501 ± 1 


21.7 ±0.1 


0.5 ±0.2 


1.16 


1.82 


1.99 



Table 2. Characteristics of microlensing candidates. 



count the experimental set up and the time sampling of 
the observations. We assume a standard model (isother- 
mal sphere with a core radius of 5 kpc) for the haloes of 
both M31 and the Milky Way. The total mass of M31 is 
assumed to be twice that of our Galaxy. MACHOs can be 
located in either haloes. Moreover, we consider also self- 
lensing due to stars in the M31 bulge or disk. We fix the 
lens masses at different values for MACHOs in the halo 
and stars in the bulge or disk. The model of the bulge is 
taken in Kent (1989), the luminosity function in Han et al. 
(1998). The luminosity function of the disk is determined 
considering two models: the one developed in Devriendt 
et al. (1999) and the model obtained considering the data 
of the solar neighborhood taken in Allen (1973) corrected 
for high luminosities (Hodge et al 1988). The results we 
obtain are almost insensitive to the particular choice be- 
tween these two models. 

We choose the mass of the lenses in the bulge to 
be 'tti'macho ~ OAMq, and the mass for MACHOs in 
the haloes equal either to m^f^QUQ — 0.5 or to 
m'lfXcHo = 0.01 Af0. In both cases about 90% of the 
expected events are due to lenses located in M31. Taking 
into account our selection criteria, the results for the ex- 
pected number of events for a halo fully composed of MA- 
CHOs, including a contribution due to lensing by stars of 
the bulge and disk of M31 of ~ 1 event independent of 
halo parameters, is ~ 4 and ~ 9, respectively. The Monte 
Carlo simulation does not yet include the effect of sec- 
ondary bumps due to artifacts of the image processing 
(alignment, seeing stabilization) and to underlying vari- 
able objects. From the data, we estimate that these effects 
reduce the number of observed events by at most 30%, and 
this for the shortest events. 

We expect, and this is confirmed by simulations, that 
the sources of most detectable events are red giants and 
have very large radii. Therefore, finite size and limb dark- 



ening effects are important, in particular for low mass 
lenses. These effects are included in the simulations^. How- 
ever, in both the real and simulated analysis, we do not 
include finite size and limb darkening effects in the ampli- 
fication light curve fitted to the events. This results in a 
loss of detection efficiency smaller than 5%, which is taken 
into account in the simulations. 

Locating our candidates in the parameter space pre- 
dicted by the simulation is more meaningful. We report 
the expected values of ii/2 and on the R magnitude at 
maximum. In Figs. |9| and |l^ we give plots of their func- 
tional relation and of their projected distributions. On 
these same plots we give the position of our selected light 
curves in this parameter space. 

Looking at the distributions we notice that for the 
^MACHO = 0.5 Mq and m'lf^cHO = 0.01 Mq cases, 80% 
of the light curves are expected to have a time width at 
half maximum ti/2 < 24 and ti/2 < 10 days, respectively 
(of course, shorter events are expected when the MACHO 
mass is smaller). In both cases, ~ 80% of the events are 
predicted with a magnitude at maximum Rmax < 21.7. 
Our candidates have ^2/2 ^ 15 days and Rmax > 21.7 
(somewhat at the limit of the expected distributions) and 
therefore most of them are probably not microlensing 
events. Still, we expect ~ 1 self-lensing event and it is 
possible that one or two of them are true microlensing 
events. In any case, from the ti/2 distribution, we are led 
to exclude that the microlensing candidate events are due 
to MACHOs of very small mass (only ^ 10% of events 
with m\f^(~,^Q = 0.01 Mq are expected to have ti/2 > 15 
days). This is indeed in agreement with the results found 
by the MACHO and EROS collaborations: they find lens 



^ The finite size effect is computed exactly from the formulae 
by Witt & Mao (1994). For the limb darkening we use an ap- 
proximation method based on Gould (1994). 
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Fig. 9. Results of Monte Carlo simulations, Ti^^^^^Q = 10^^ Mq case. The scale on the y coordinate of the two 
distributions on the right are in arbitrary units. The dots (left) give the position of the candidates as labelled in Table 
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Fig. 10. Results of Monte Carlo simulations, ra^M^^CHO ~ ^.5 Af© case. The scale on the y coordinate of the two 
distributions on the right are in arbitrary units. The dots (left) give the position of the candidates as labelled in Table 



masses in the halo within the range 0.2-0.6 Mq (Alcock 
et al. 2000; Lasserre et al. 2000). 

We are not yet in a position to tell what kind of vary- 
ing objects generate our events if they are not due to mi- 
crolensing. They may be irregular or long period variable 
giants, but only a longer time baseline, and/or observa- 
tions of the object at minimum light, will allow us to con- 
clude. 

The MDM analysis is not yet complete. The results 
from the analysis of data acquired in the other field (lo- 
cated on the opposite side with respect to the major axis 
of M31 of the Target field analysed here) and results from 
new observations scheduled for the fall 2001 will give us 



the opportunity to gain further insight into the still open 
question of the composition of dark haloes. At the present 
time, with the caution suggested by the just mentioned 
problems, the analysis discussed in this paper tends to 
confirm that only a minor fraction of dark matter in the 
galactic haloes is in form of MACHOs within the mass 
range 0.01-0.5 M© under the assumption of a standard 
model of the halo and given source luminosity functions. 
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